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Abstract 
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renormalization-group flow and the phase transitions of random tensor networks. 
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1 Introduction 


Wilson’s renormalization group is an essential and pedagogical tool in modern theoretical 
physics. Once a renormalization-group flow in a parameter space is given, one can read off 
relevant degrees of freedom at each step of coarse graining through change of parameters, and 
understand the phase structure in principle. Therefore, a renormalization-group flow gives 
us a quantitative and qualitative picture of a system concerned. The aim of this paper is 
to define a renormalization procedure and derive the corresponding flow equation for random 
tensor networks, in particular for those proposed as Feynman-graph expressions M, through 
the use of the canonical tensor model (CTM, for short). 

First of all, CTM has been introduced by one of the authors as a model of quantum gravity 
by considering space-time as a dynamical fuzzy space [SUSIE]- CTM is a tensor model in the 
canonical formalism, which has a canonical conjugate pair of rank-three tensors, M a b c , P a b c 
( a,b,c = 1, 2, • • • , TV), as dynamical variables. This interpretation of tensorial variables in 
terms of a fuzzy space is different from the one made by original tensor models. Historically, 
tensor models have been introduced as models of simplicial quantum gravity in dimensions 
higher than two [SI El ED]; although the original tensor models have some drawbacks, tensor 
models as simplicial quantum gravity are currently in progress as colored tensor models mm 
(See [13]-[53] for recent developments.). In CTM, N, the cardinality of the rank-three tensors, 
may be interpreted as the number of “points” forming a space, while physical properties 
of space-time such as dimensions, locality, etc. must emerge from the collective dynamics 
of these “points.” So far, the physics of the small- TV CTM is relatively well understood: the 
classical dynamics of the N — 1 CTM agrees with the minisuperspace approximation of general 
relativity in arbitrary dimensions [25]; the exact physical states have been obtained for N = 2 
in the full theory [26] m\ and for N = 3 in an S' 3 -symmetric subsector [27]; intriguingly, 
physical-state wavefunctions, at least for N = 2,3, have singularities where symmetries are 
enhanced [27]. However, similar brute-force analysis as above for N > 3 seems technically 
difficult because of the huge number of degrees of freedom of the tensorial variables, although 
in order to capture, for instance, emergence of space-time from CTM, the large-iV dynamics 
is supposed to be important. Thus, for the purpose of handling large- N behaviors of CTM, 
the present authors have proposed the conjecture that statistical systems on random networks 
[DUSES], or random tensor networks, are intimately related to CTM [3]: the phase structure of 
random tensor networks is equivalent to what is implied by considering the Hamiltonian vector 
flow of CTM as the renormalization-group flow of random tensor networks. This conjecture 
has been checked qualitatively for N = 2 [3]. I 11 fact, as more or less desired, random tensor 
networks turn out to be useful to find physical states of CTM with arbitrary N: some series of 
exact physical states for arbitrary N have been found as integral expressions based on random 
tensor networks [27]. 

In this paper, we prove the fundamental aspect of the above conjecture: we show that the 
Hamiltonian vector flow of CTM can be regarded as a renormalization-group flow of random 
tensor networks for general N. Here the key ingredient is that the Lagrange multipliers of the 
Hamiltonian vector flow are determined by the dynamics of random tensor networks in the 
manner given in this paper. This is in contrast with the previous treatment for N — 2, in 
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which the Lagrange multipliers are given by a “reasonable” choice [3j. In fact, the previous 
treatment turns out to have some problems for general N, as being discussed in this paper. 

This paper is organized as follows. In section [21 we review CTM and random tensor net¬ 
works. We argue our previous proposal [3j on the relation between CTM and random tensor 
networks, and its potential problems. In Section [3j we propose a renormalization procedure 
for random tensor networks based on CTM, and derive the corresponding renormalization- 
group flow. In Section we compare our new and previous proposals with the actual phase 
structures of random tensor networks for N = 2,3. We find that the new proposal is con¬ 
sistent with the phase structures, while the previous one is not. In Section [5j we discuss the 
asymptotic behavior of the renormalization-group flow, and clarify the physical meaning of 
the renormalization parameter. In Section [6l we provide a general theorem which relates dis¬ 
continuity of the renormalization-group flow and the phase transitions of the random tensor 
networks. Section [7] is devoted to summary and discussions. 


2 Previous proposal and its problems 


In this paper, we consider a statistical system [310] parameterized by a real symmetric three- 
index tensor P abc (a,b,c — 1,2,..., iV) J3 Its partition function is defined b}0 


Z n {P) = J d(t> (P0 3 ) n e~ n ^ 2 , 


( 1 ) 


where we have used the following short-hand notations, 

/ N ,00 

d(j) = I d(j) a i P(j) = Pabc^a^b^ct 4* * = (fra&a- 

o=l 


( 2 ) 


By using the Wick theorem, the Gaussian integration of </> in ([TJ) can be evaluated by the 
summation over the pairwise contractions of all the 0’s in (P0 3 ) n . Then the partition function 
(JT|) can graphically be represented by the summation over all the possible closed networks of 
n trivalent vertices. In each of such networks, every vertex is weighted by P a bc , and the 
indices are contracted according to the connection of the network, as in Fig{T] Therefore, 
since the pairwise contractions are taken over all the possible ways, the statistical system 
represented by (JTJ) can be regarded as random tensor networks of n trivalent vertices. In 
general, such a network may contain disconnected sub-networks, but this probability vanishes 
in the thermodynamic limit n —)■ ooji] 

*In this paper, the tensor variable of the statistical system is denoted by P for later convenience, instead 
of M used in the previous papers E |4j. 

tFor later convenience, the normalizations of the partition function and cf> are taken differently from those 
in the previous papers EE- This does not change the physical properties of the statistical system. 

* This graphical property can be checked by considering a sort of “grand” partition function Z(P) given by 
a formal sum of JT]) over n with n-dependent weights, and comparing explicitly the perturbative expansions in 
P of Z{P ) and those of log Z(P) for N = 1. The former corresponds to the sums of the networks which may 
contain disconnected sub-networks, while the latter to connected networks only. See {3] for details. 
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Figure 1: 



For example in N = 2, (pQ) gives the partition function of the Ising model on random 
networks M, if one takes 

2 

Pabc = ^2 R aiRbiRcie H(Ti , (3) 

1=1 


where a l represents the spin degrees of freedom taking a\ = 1, <72 = — 1, H is a magnetic held, 
and R is a two-by-two matrix satisfying 

(R T R) ij = e J ^ : (4) 


with J giving the nearest neighbor coupling of the Ising model. For a ferromagnetic case 
J > 0, there exists a real matrix R satisfying (HI) . 

The partition function (JT() is obviously invariant under the orthogonal transformation L G 
O(N), which acts on P as 

Pabc ^-‘ aa 'I J bb'I J cc'P a 'b’c'i (b) 

since the transformation can be absorbed by the redefinition <// = (p a 'L a i a . In addition, the 
overall scale transformation of P, 

( 6 ) 


with an arbitrary real number -0, does not change the properties of the statistical system, 
since this merely changes the overall factor of dU). For example, for N = 2, these invariances 
allow one to consider a gauge, 


Pill — 1, P 112 — 0, P\22 ~ Xl, P222 ~ X2, 


(7) 


with real x % . 

The free energy per vertex in the thermodynamic limit can be defined by 


f(P) 


lim 

n —>00 


log Z 2n (P) 


2 n 


J 


( 8 ) 


where we have considered only even numbers of vertices, since an odd number of trivalent 
vertices cannot form a closed network. The phase structure of the statistical system can be 
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Figure 2: The horizontal and vertical axes represent X\ and x 2 of the gauge (J7J), respectively. 
The solid lines describe the phase transition lines of the random tensor networks with N = 2, 
which can exactly be obtained in the thermodynamic limit [3]. The arrows describe the 
Hamilton vector flow (12TD of CTM with J\f a = P a bb- 



investigated by studying the behavior of the free energy (jSJ). For the case of N = 2, the phase 
transition lines of the free energy ([HD in the gauge (ED are shown by solid lines in Fig|2] [3] 0]. 
The transitions are first order, except at the Curie point, {xi,x 2 ) = (0.5,0), where the first 
derivatives of /(P) are continuous, but the second ones are not p.:, HI 28j. In fact, for arbitrary 
N, the free energy can exactly be obtained by applying the Laplace method to evaluate ((HD 
011 ]. The result i^§ 

/(P)=Min*/(P,0) = /(P,$, (9) 


where 


f{PA) = <i > 2 


1 

2 


log 


3p 


(py) 


( 10 ) 


and (ft is defined so as to minimize /(P, (ft) as a function of (ft for given P, namely, (ft is one of 
the solutions to the stationary condition, 


<9/(P» 

d(ft a 


2 (fta - 


3Pabc(ftb(ftc 

P03 


0. 


( 11 ) 


Tf P is symmetric under part of the O(N) transformation (0, the perturbations of the integrand of 0 
around (ft = (ft contain some zero modes, and the application of the Laplace method will require extra treatment 
to integrate over the symmetric directions. However, this integration is obviously finite, and will only generate 
corrections of the free energy higher in 1/n, which do not affect the thermodynamic limit. Thus, the exact 
free energy in the thermodynamic limit is given by 0 for the whole region of P including symmetric points. 
Note that this argument may change, if one takes a double scaling limit accompanied with iV —> oo, which, 
however, is not considered in this paper. 
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The most important implication of the previous paper [3] was that the phase structure of 
the Ising model on random networks (more exactly, random tensor networks with N — 2) in 
Fig[2] can be derived from the Hamilton vector flow of CTM for N — 2, if one regards the 
Hamilton vector flow as a renormalization-group flow, as shown in Figj2j This is surprising, 
since CTM was proposed aiming for quantum gravity, and there exist no apparent reasons for 
CTM to be related to statistical systems on random networks. CTM is a totally constrained 
system with a number of first-class constraints forming an algebra which resembles the Dirac 
algebra of the ADM formalism of general relativity [29]. In the classical case, the constraints 
are given by 


U a 2 P ah cPbdeM cde , ( 12 ) 

J o,6] = - (PacdMbcd — PbcdM ac d ) , (13) 

P = X PabcM a bc , (14) 

O 

where J and V are the kinematical symmetry generators corresponding to the SO(N ) (J5]) 
and the scale (J6j) transformations, respectively, and and J may be called Hamiltonian and 
momentum constraints, respectively, in analogy with general relativity [29]. Here the bracket 
for the indices of J symbolically represents the antisymmetry, J [ ab ] = —J[b a ], and M is the 
canonical conjugate variable to P defined by 

{■M a bci Pdef} g ^ ^ Aiafr/) ^ha(e) ^co(f) • {-M a bci -Mdef} {Pabci Pdef} 0, (1^) 

a 

where { , } denotes the Poisson bracket, and the summation over a runs over all the per¬ 
mutations of d, e, f to incorporate the symmetric property of the three-index tensors. The 
constraints form a first class constraint algebra, 


{WK 1 ),WK 2 )} = ^([« 1 .f 2 ]), 


{JM, m)} = \u (,}(): 

(16) 

{j(y),j(K 2 )} = lj(bV i), 

(17) 

{Z>, «(«)} = 

(18) 

= o, 

(19) 


where U(£) = J(v) = V[ab]J[ab] and £ ab = P a bdc- 

The Hamiltonian of CTM is given by an arbitrary linear combination of the constraints as 


H — J\f a P-a + A f[ab]<J[ab\ + J\fP, 


( 20 ) 
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where A/"’s are the multipliers, which may depend on P in the context of this paper, considering 
a flow in the configuration space of P . Then, the Hamiltonian vector flow is given by0 

^Pabc={H,P abc }, (21) 

where s is a fictitious parameter along the flow. In the previous paper [3], which compares 
CTM with the random tensor networks for N = 2, the multiplier A f a is chosen to be 

K = Pabb , (22) 

based on that this is the simplest covariant choice. The other multipliers Af[ ab ] and A/", related 
to the symmetry generators, are chosen so that the Hamilton vector flow fl2T|) keeps the gauge 
condition (|7|). Indeed the flow in Fig J2] has been drawn with these choices. One can also check 
that other covariant choices such as J\f a = P a b c PbdePcde do not change the qualitative nature 
of the flow and therefore the coincidence between the phase structure of the random tensor 
networks and the one implied by CTM with N = 2. 

Though the coincidence is remarkable for N = 2, from further study generalizing gauge 
conditions and values of N, we have noticed that there exist some problems in insisting the 
coincidence as follows: 

• First of all, no physical reasons have been given for the coincidence. A primary expec¬ 
tation is that there exists a renormalization group procedure for statistical systems on 
random networks, and the procedure is described by the Hamiltonian of CTM in some 
manner. However, it is unclear how one can define a renormalization group procedure for 
statistical systems on random networks, which do not have regular lattice-like structures. 

• As will explicitly be shown later, in the case of N = 3, the phase transition lines deviate 
from the expectation from the Hamilton vector flow of CTM. What is worse is that 
different choices of A f aj such as P abb or PabcPbdePcde, give qualitatively different Hamilton 
vector flows, which ruins the predictability of the transition lines from the flow. 

• In Figj2]for N = 2, on the phase transition lines, the flow goes along them, and there 
exist a few fixed points of the flow on the transition lines. The fixed point at (0,1) is a co¬ 
dimension two fixed point, and the associated phase transition is expected to be of second 
order rather than first order, if the flow is rigidly interpreted as a renormalization-group 
flow and we follow the standard criterion [SU]. This is in contradiction to the actual 
order of the phase transition. This contradiction is more apparent in the diagram in 
another gauge in Section [0 

The purpose of the present paper is to solve all the above problems, and to show that 
CTM actually gives an exact correspondence to random tensor networks. It turns out that M a 
should not be given by any “reasonable” choices as above, but should rather be determined 
dynamically as N a ~ {4>a) to be discussed in the following sections. Then, we can show that 
the Hamiltonian of CTM actually describes a coarse-graining procedure of random tensor 
networks, and that the Hamilton vector flow is in perfect agreement with the phase structure 
irrespective of values of N. 

^The direction of the flow is chosen in the manner convenient for later discussions. 
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3 Renormalization procedure and renormalization-group 
flow 

In this subsection, we discuss a renormalization group procedure of the random tensor net¬ 
works, and obtain the corresponding renormalization group flow. 

Let us consider an operator O which applies on Z n (P) as 



(23) 


By using (fl2li and (fT5lh and performing partial integrations with respect to </>, one can derive 



3n+3+iV 



(24) 


where (4> a )n is an expectation value defined by 



(25) 


and the numerical factor in the first term of (1241) is due to the rescaling of (j) for ruj) 2 —> ( n+l)q i 2 


in the exponential. 

Here (HM]) and (125|) must be used with a caution. If taken literally, since Z n=odd = 0, (OH) 
and (l25j) do not seem useful by themselves. The reason for Z n=odd = 0 is that the contributions 
at (j) = ±u with arbitrary v cancel with each other in the integration of ([[]). To avoid this 
cancellation and make (1241) and (1251) useful, let us consider a finite small region in the space 
of (j) around one of the solutions </> which minimize m- For later convenience, we take the 
sign of cj) so as to satisfy 


P(f) 3 > 0. 


(26) 
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This can be taken, because, if P0 3 = 0, /(P, 0) in (TTOT) diverges and cannot be the minimum, 
and f(P, 0) = f(P, —0)- Especially, should not contain the other minimum 0 = — 0 . Then, 
let us consider a replacement, 


Z n {P )-»• / d0 (P0 3 )"e-^ 2 . (27) 

In the thermodynamic limit n —> oo, the integral (12711 is dominated by the region with width 
A0 ~ l/\/n around 0 = 41- Therefore, the expression (T27|) approaches e~ n PP in the ther¬ 
modynamic limit, irrespective of even or odd n. Moreover, since the integrand of (1271) damps 
exponentially in n on the boundary of r^, the corrections generated by the partial integrations 
carried out in the derivation of (12211 are exponentially small. Thus, (122)1 is valid up to expo¬ 
nentially small corrections in n after the replacement (1271) . Thus, for n large enough, we can 
safely use (124)1 and (125)1 as if they are meaningful irrespective of even or odd n. 

Taking into account the discussions above, we can put (0 a ) n — > 0 a and Z n —> e~ n ^ P) in 
(122)1 for n 1. Therefore, in (122)1 . the first term dominates over the second term, and one can 
safely regard O as an operator which increases the size n of networks. To regard this operation 
as a flow in the space of P rather than a discrete step of increasing n, let us introduce the 
following operator with a continuous parameter s, 

R(s ) = e s °. (28) 

If we consider n 1, we can well approximate the operation O with the first term of (12211 as 
explained above, and one obtains 


qTn 

R(s)Z n (P) = V — -O m Z n (P) 
ml 

m= 0 

^ s m (n + m- l)! c _ (n+m)/(P)+ 3m. 

Z-^ 3 m m\(n — 1)! 

m =0 v ’ 

By increasing s, the right-hand side is dominated by larger networks, and diverges at s = Soo, 
which is given by 



Soo = 3 exp (f(P) - 0 . (30) 

On the other hand, in the thermodynamic limit, the left-hand side of (129)1 can be computed 
in a different manner. In the thermodynamic limit, 0 a can be replaced with the mean value 
0 a , and the operator O can be identified with a first-order partial differential operator, 

O —> Or) = 0 o P a = ~([)aPabcPbdeD P de: (31) 

Hlf P is symmetric under part of the SO{N ) transformation ([5]), an extra care will be needed as discussed 
in a previous footnote. However, this does not change the following argument in the thermodynaic limit. 




where D^ bc is a partial derivative with respect to P abc with a normalization, 


DabcPdef — { M abc , Pdef} 


1 
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^ ^ ^acr(d) $bo(e) ^ccr(f) • 


( 32 ) 


Here note that <p is a function of P determined through the minimization of (fTOlh Then, the 
expression in the left-hand side of (1291) is obviously a solution to a first-order partial differential 
equation, 


“ ° D ) R ^ Zn ^ = °- ( 33 ) 
The solution to ()33l) can be obtained by the method of characteristics and is given by 

R(s)Z n (P) = Z n (P(s )), (34) 

where P(s) is a solution to a flow equation, 

-J -Pabc(s) = 0 D Pabc(s ) 

as 

g fidPdae(s') Pebc^S ) T (pdPdbe(s') Peca^S ) T (p dPdcei^S) -fea&(®)) j (35) 

Pabc(s = 0) = P a6c , (36) 

where <p must be regarded as a function of P(s) by substituting P with P(s). 

Here we summarize what we have obtained from the above discussions. R(s)Z n (P) can 
be evaluated in two different ways. One is (|29]1 . a summation of random tensor networks, the 
dominant size of which increases as s increases, while P is unchanged. The other is (l34li , where 
P(s) changes with the flow equation (j35]h while the size of random networks is unchanged. 
This means that the change of the size of networks can be translated into the change of P. 
Therefore the flow of P(s) in (1351) can be regarded as a renormalization-group flow of the 
random tensor networks, where increasing s corresponds to the infrared direction. 

The above derivation of the renormalization group flow uses the particular form of R 
in (H2D- Since, in general, there exist various schemes for renormalization procedures for 
statistical systems, one would suspect that there would be other possible forms of % which 
describe renormalization procedures for random tensor networks. However, this is unlikely, 
and the form (fl2|) would be the unique and the simplest. The reason for the uniqueness is 
that, as outlined in Section [2l the algebraic consistency of R with the O(N) symmetry, which 
is actually the symmetry of random tensor networks in the form (jT]) , requires the unique form 
(J2D under some physically reasonable assumptions |6j. On the other hand, the reason for the 
simplest can be found by considering the diagrammatic meaning of the operation R in ([23]) . 
R acts on a vertex as 

{( p a R a ,P(p 3 } = ^(p a (pbPabcPcde<Pd<Pe , (37) 
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Figure 3: The diagrammatic representation of the operation {0 a "H a ,P0 3 }. 


and hence can be regarded as an operator which inserts a vertex on an arbitrary connection 
in a network (See Figj3|). This is obviously the most fundamental operation which increases 
the number of vertices of a network. 

Here we comment on our new proposal in comparison with the previous one |!3]. Our main 
claim is that the multiplier should take J\T a = 4>a rather than “reasonable” choices such as 
J\f a = P a bb , PabcPbedPcde, etc., taken in the previous proposal. With J\f a = 0 a , the Hamiltonian 
vector flow is uniquely determined by the dynamics, while “reasonable” choices are ambiguous. 
Even if ambiguous, there are no problems in the N — 2 case, since there are no qualitative 
changes of the flow among “reasonable” choices, and the phase structure can uniquely be 
determined from the flow. However, as will be shown in Section SI this is not true in general 
for N > 2. In fact, N = 2 is special for the following reasons. It is true that (f) a can well be 
approximated by ~ const.P a bb near the absorbing fixed points in Fig|2j because all of them 
can be shown to be gauge-equivalent to Pm = 1, others = 0. This means that at least an 
approximate phase structure can be obtained even by putting Af a = P a bb ~ 4>a- In addition, 
what makes the N — 2 case very special is that the phase transition lines are the fixed points 
of the Z 2 symmetry corresponding to reversing the sign of the magnetic held of the Ising model 
on random networks. Therefore, the phase transition lines are protected by the symmetry, 
which stabilizes the qualitative properties of the how under any changes of Af a respecting the 
symmetry. 

Finally, we comment on an equation which can be derived from (170) in the thermodynamic 
limit. By putting Z n (P) ~ to (EMI) , one can derive 

4> a P abc P bde D^ de f{P) = - 2 -e~l^ +3 Z (38) 

In fact, one can directly prove (138]h By using (J§D and (HUD , the left-hand side of (155jl is given 
by 




Ip p (n p I 

(PaPabcPbde ( PcdeVg) o 7 

d( Pg 

2 

—3 p *’ 


0 a Pabc Pbde 



P(J) 3 


(39) 


where we have used (fill) . This coincides with the right-hand side of (|38|) . because of the choice 
(T 26 |) and (j) 2 = |, which can be obtained by contracting (flTTl with 0 a . 
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Figure 4: The phase transition lines of the random tensor network with N — 2 in the gauge (14(1 
are shown as solid lines. The horizontal and vertical axes indicate aq and x 2 , respectively. The 
phase transitions are first-order except for the endpoint of a line located around ( Xi,X 2 ) = 
(0.2,1.2). This is the Curie point, which is gauge equivalent to (aq,^) = (0.5,0) in the 
gauge ([7]). The left figure describes the Hamilton vector flow based on our former proposal, 
Af a = P a bbi while the right figure describes it based on our new proposal, Af a = 0 a . A locus of 
gauge singularities is located at x\ = 0.635. Gauge singularities are not physical, and the free 
energy has no singular behaviors there. Some details are given in Appendix IA.1I 


4 Comparison 

In this section, we will check the proposal of this paper in the cases of N = 2, 3 by comparing 
with the phase structures of the random tensor networks. 

Let us first consider the N = 2 case with a gauge, 

7*111 = 1) 7*112 = 0.3, 7*122 = x li 7*222 = x 2 j (40) 

as a typical example. The difference from (|7j) is the gauge fixing value of Pm- One can 
obtain the phase structure in the parameter space of (aq,^) by studying the free energy ([9]). 
Alternatively, one can apply the 0(2) and scale transformations, (J5]) and |6|, on F so that the 
phase structure in the gauge ([7]) given in Figj2] is transformed to that in the gauge (1401) . In 
either way, one can determine the phase structure in the new gauge, and the result is given in 
FiglU Here we draw the Hamilton vector flows for Af a = P a bb, based on the former proposal, 
and Af a = 4> a , based on our new proposal, in the left and right figures, respectively. 

The rough features of the two flows based on the different proposals seem consistent with 
the phase structure: the flows depart from the transition lines, and go into the same absorption 
fixed points. This was the main argument in our previous paper [3]. However, there are some 
physically important differences in details between the left and the right figures. In the left 
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figure, on the phase transition lines, the flow is going along them. Moreover, there exist a few 
fixed points of the flow on the transition lines at {x\,X 2 ) 0.2,1), (-0.6, 0.5), (-0.4,-1) 

in the left figure. If the flow is strictly interpreted as a renormalization-group flow, the phase 
transition line on the righthand side of the fixed point near (— 0 . 2 , 1 ) is expected to be of 
second order, rather than first order, since the points on the both sides of the transition line 
in its vicinity flow to the same fixed point near (0.1,0) without any discontinuities. On the 
other hand, in the right figure, the flow has discontinuity on the transition lines, except for 
the Curie point at the endpoint of the transition line. Thus, the flow based on our former 
proposal clearly contradicts the actual order of the phase transitions, while the one based our 
new proposal is in agreement with it, i.e. first order except for the Curie point. 

An interesting property of the flow is that it does not vanish even on the Curie point, 
as can be seen in the right figure of FigJT] and can also be checked numerically. This seems 
curious, because the second derivatives of the free energy contain divergences on the point. In 
a statistical system on a regular lattice, such divergences originate from an infinite correlation 
length. Therefore, such a point will typically become a fixed point of a renormalization-group 
flow. On the other hand, the correlation length of the Ising model on random networks is 
known to be finite even on the Curie point [25]. This means that, even starting from the Curie 
point, a renormalization process will bring the system to one with a vanishing correlation 
length. This implies that the Curie point cannot be a fixed point of a renormalization-group 
flow, and this is correctly reflected in the fact that our flow does not vanish on the Curie point. 

Let us next consider the N = 3 case. There seem to exist too many parameters to treat 
this case in full generality. So let us specifically consider a subspace parametrized by 

Piii = 1, Pijj = x i, P123 = x 2, (i j), (41) 

which is invariant under the S 3 transformation permuting the index labels, 1,2,3. Through 
numerical study of the free energy (|U]) (and some analytic considerations), one can obtain the 
phase structure shown in Fig]5] I 11 the indicated parameter region, there exist two regions of 
an S 3 symmetric phase labeled by S with 0i = 0 2 = 03- There also exist two distinct non- 
symmetric phases labeled by NS1 and NS2. At any point in the two regions, the minimization 
of the free energy ([TUI) has three distinct solutions of non-symmetric values, </>i 7 ^ 0 2 = 03, 
02 7 ^ 03 = 03 7 ^ 0i = 02 , and hence three distinct phases coexist in these regions. When 

S 3 symmetric subspace (ITT]) is extended to more general cases, each of NS1 and NS 2 becomes 
the common phase boundary of the three phases. 

The flow in the left figure of Fig]5] is drawn based on our previous proposal J\f a = P a bb- 
There, the flow is not in good agreement with the phase structure, though it seems to capture 
some rough features. We tried other possibilities such as J\f a = PabcPbdePcde, etc., but the 
flow depended on the choices, and no good agreement could be found. On the other hand, 
in the right figure based on our new proposal A f a = 0 a , the flow in the symmetric region S is 
consistent with the phase structure: the flow departs from the transition lines, and, since it 
does not vanish on the lines, the order is expected to be first order. This is in agreement with 
the property of the free energy except for a point (xi,x 2 ) = (1/4,— 1 / 8 ). At this point, the 
free energy is continuous in the first derivatives, but singular in the second. However, since 
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Figure 5: The phase transition lines of the random tensor network with N — 3 in the subspace 
(14B are shown as solid lines. The horizontal and vertical axes indicate X\ and X 2 , respectively. 
There exist four regions separated by the phase transition lines in the indicated parameter 
region. The phases labeled by S have S 3 symmetric expectation values </>i = </>2 = 4 > 3 - There 
are two other regions labeled by NS1 and NS2, which have non-symmetric expectation values. 
The phase transitions are first-order except for the meeting point of S, NS1, NS2 at ( xi,x 2 ) = 
(1/4, —1/8). On this point, the free energy is regular in the first derivatives, but singular in 
the second. On the phase transition line between NS1 and NS2, P is symmetric under an 
SO (2) transfomation [27] . In the left figure, we draw the Hamilton vector flow based on our 
former proposal Af a = P a bb, while the right one is based on our new proposal Af a = (pa- 
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the flow does not vanish on the point, the correlation length is expected to be finite. This is 
similar to the Curie point of the N — 2 case. 

In the non-symmetric phases, NS1 and NS2, the expectation values are not S 3 symmetric. 
Therefore, the flow has generally directions away from the S 3 symmetric subspace, and cannot 
be drawn on the figure. To also check the consistency of the flow in these regions, it would be 
necessary to extend the parameter region. This would require to take a different systematic 
strategy for consistency check to avoid too many parameters. In Section [HI we will take another 
way of consistency check by proving a theorem relating the renormalization group flow and 
the phase transitions of the random tensor network. 


5 Asymptotic behavior of the flow 

In Section [3j we argued that O = (paHa provides a renormalization procedure for the random 
tensor network. As can be seen from (1291) and (1341) . P(s) diverges in the limit s —» Soo in (1301) . 
On the other hand, in the numerical analysis of Section [4J, P is kept normalized as (1401) and 
m by appropriately tuning the multiplier A f for the scale transformation V (with J\f a bJab as 
well). As in FigQ] and G3 one can find fixed points of the flows in the limit s —>- 00, where s 
denotes the fictitious parameter parameterizing the normalized flows. In this section, we will 
show that these two limits of s and s are physically equivalent. 

Let us first show the divergence in s —> more directly. Since (p 2 = § from dll]) , 

log P(p 3 = —f(P) + f. Then, by using (13TT) and (139|) . one obtains 

^-log Pfi= l -Pj>\ (42) 

as 6 

where P is meant to be P(s), and hence Pep 3 is regarded as a function of s. The solution to 

(USD is 


p (p 3 


1 

1 ~ ’ 

p **u 5 


(43) 


which indeed diverges at s = ^4— = s^. Since (p is normalized by (p 2 = |, the divergence 

^ I 3=0 

of (1T5D can be translated to the divergence of P{s) with a behavior, 


P{*)~ 


const. 



(44) 


or higher order in the case that some components of (p vanish in s —> Sqo- 

Now let us compare the two flows, unnormalized and normalized ones. For notational 
simplicity, let us denote the three indices of P a b c by one index as P % . The flow equations in s 
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and s can respectively be expressed as 


^-P,(s) = i a (P(s)) g a ,(P(s)). 
ds 

(45) 

I 'zPi(§ ) = MP(s))gai(P(s)) -AfPiCs), 

(46) 

o 

> 

'Zu 

o 

(47) 


where A f generally depends on P t (s), and g a i(P ) are the quadratic polynomial functions of P, 
which can be read from (135D . The last term of the second equation comes from J\fV in (120]) , 
and is assumed to be tuned to satisfy a gauge condition normalizing P(s). Here we ignore the 
SO(N ) generators, J7ab, for simplicity, but it is not difficult to extend the following proof to 
include them. 

The physical properties of the random tensor network do not depend on the overall scale 
of P. So let us define the relative values of P(s) and P(s) as 


Qi(s) 


Pj(s) 

Po(sY 


Qi{s ) 


A(5) 

PoCsY 


(48) 


where Po(s) (Po(s)) is taken from one of Pi(s) (resp.Pj(s)), or a linear combination of them. 
From (HID , it is obvious that 0(P) and 0(P) actually depend only on Q and Q , respectively. 
Then, from Ql5lh 

Id 1 d Pi 

Po ds 1 Po ds Po 

_ 4>a{Q) 9ai(P) _ 4>a(Q) 9ao(P)Pi 

p2 p3 

= 4>a(Q) (g a i(Q) - 9ao(Q)Qi) ■ (49) 

In the same manner, 

= MQ) ( 9ai{Q) - 9ao(Q)Qi ) • (50) 

Note that the last term of (14611 does not contribute to the flow equation of Q. Since the initial 
condition (HTD implies Qi(s — 0) = Qi{s — 0), and the righthand sides of (H5]l and (HUD are 
identical, the flow equations, (H9j) and (150|1 . describe an identical flow with a transformation 
between the fictitious parameters, 


Po(s)ds = Po(s)ds. 


(51) 


As discussed above, the typical behavior of Po(s) is (144]) . while P 0 (s) is assumed to remain 
finite near an absorption fixed point. In such a case, (l5T]l implies 

s ~ —const. log(soo — s). (52) 
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Therefore, the limits of s —y s^ and s —> oo are physically equivalent. As can be checked easily, 
this physical implication does not change, even if we consider the case that Pq(s) diverges with 
an order higher than (144 ft . 

To investigate the physical meaning of the fictitious parameter s, let us estimate (l29]i near 
S rs ~' Soo. We obtain 


s m (n + m- 1)! c -( n+m )/ (P ) + 3m 
^ 3 m m\(n — 1)! 

m=0 v 7 


e -nf(P) “ 

7 -ttt > (m + 1) (m + 2) • • • (n + m 

(n — 1 ! 

^ 7 m— n 


= const.- 


d 


m =0 
n— 1 00 




m=0 


"- 1 ^ ( Soo 


m+n—1 


= const.- 


d n-l 8 n-l 


ds 


n— 1 


^oo ^ 


const, (soo — 


1) 



m 


(53) 


Then the average size of networks can be estimated as 

(n + m) ~ n + s-f- logfsoo — s) -n ~ nS °° . (54) 

as Soo — s 

Therefore 


s ~ const, log (Average size). (55) 

This means that s corresponds to the standard renormalization-group scale parameter often 
denoted by log A in field theory. 


6 Discontinuity of the renormalization-group flow and 
phase transitions 

In Section [H we see that the renormalization-group flow has discontinuity on the first-order 
phase transition lines in the examples of the random tensor networks. In this section, we will 
prove a general theorem on this aspect. 

By using the free energy in the thermodynamic limit ()9j) . the stationary condition (HIT) and 
the flow equation (j35]h we can prove the following theorem. 

Theorem: The following three statements are equivalent. 

(i) The first derivatives of f(P) are continuous at P. 

(ii) (f) is continuous at P. 

(iii) The renormalization-group flow is continuous at P. 


16 







Proof: 

Let us first prove (i) =>- (ii). From (]9]h the first derivatives of /(P) are given by 


D abJ( P ) = 



( 56 ) 


where we have neglected the contributions from the P-dependence of (ft, since (ft satisfies the 
stationary condition (fill) . By contracting a pair of indices in (1561) . one obtains, 


D P aab f(P) 



3 (ftb 
2P0 3 ’ 


(57) 


where we have used (ft 2 = | derived from CD- Here note that P(ft 3 is continuous at any P, 
because the free energy /(P) itself in (jHJ) is continuous at any F0, and also (ft 2 = |. Therefore, 
if (i) holds, (1571) is continuous and hence (ft is continuous; (i) =>- (ii) has been proven. 

The reverse, (ii) =>■ (i), is obviously true from (1561) . Therefore, the statements (i) and (ii) 
are equivalent: (i) (ii). 

Next, as for (ii) (iii), it is obvious that, if (ft is continuous, the renormalization group 
flow (j35]) is also continuous. 

Finally, let us prove (iii)=^(ii), which will complete the proof of the theorem. To prove this, 
we will show that there is a contradiction, if we assume both (iii) and that (ft has discontinuity 
on P. 

Let us suppose that there is discontinuity of (ft at a point P. Then, from the definition of 
(ft, there exist multiple distinct solutions of (ft to (ITT]) which give the same minimum of (fTOl) at 
P. Let us take any two of them, (ft + and (ft~. As shown above, P0 3 is continuous at any point, 
which means 


A = P(i+) 3 = P(i-f, 


(58) 


where the value is denoted by A for later usage. Here note A ^ 0, since, otherwise, (fTOl) 
diverges and cannot be the minimum. Then, since (ft ± both satisfy (ITT]) , we obtain 

(fta = 7^ P abc(ftb(ftt, ( 59 ) 

A a = ^ Pabc (AbA c + 2A60") , (60) 

where A = (ft + — (ft~, and (160]) has been obtained by considering the difference of the two 
equations in (|59|) . Note A ^ 0, if there exist multiplicity of the solutions. 

On the other hand, the assumption (iii) and (135|) imply 


ds 


-P 


abc 


<j>=(t>+ 


_± P 

7 1 abc 

as 


OC A dPdaePebc + A dPdbePeca + A dPdcePeab ~ 0. (61) 


**This can be proven by using that f(P) is the minimum of (1101) . which is a continuous function of (ft and P. 
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Then, by contracting (RjTj) with three 0+’s or 0 ’s, and using (I59j) . we obtain 

P abc A a 0+0+ = 0. (62) 

Finally, by contracting (I60p with A, we obtain 

A 2 = ^ P abc (A a A b A c + 2A a A b (p~) 

= ^PabcAa (0+0+ - 0^0") 

= 0, (63) 

where we have used (l62]h This concludes A = 0, which contradicts the initial assumption of 
the existence of discontinuity of 0 . Consequently, we have proven the equivalence of (i), (ii) 
and (iii). I 


By taking contrapositions, a corollary of the theorem is given by 
Corollary 1: The following three statements are equivalent. 

(i) P is a first-order phase transition point. (Not all of the first derivatives of f(P) are 
continuous.) 

(ii) 0 is not continuous at P. 

(iii) The renormalization-group flow is not continuous at P. 

Another corollary of physical interest is 

Corollary 2: If P is a phase transition point higher than first-order, the renormalization- 
group flow is continuous at the critical point. 

The qualitative behavior of the N = 2 renormalization-group flow shown in the right fig¬ 
ure of Figfl] respects the theorem and corollaries as it should be: Corollary 1 is realized on the 
phase transition lines, and Corollary 2 on the Curie point. 


7 Summary and discussions 

In the previous paper [3], it has been found that the phase structure of the Ising model on 
random networks (or random tensor networks with N = 2) can be derived from the canon¬ 
ical tensor model (CTM), if the Hamilton vector flow of the N = 2 CTM is regarded as a 
renormalization-group flow of the Ising model on random networks. This was a surprise, since 
CTM had been developed aiming for a model of quantum gravity in the Hamilton formalism 
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[3 IS, 5]- Considering the serious lack of real experiments on quantum gravity, the aspect that 
CTM may link quantum gravity and concrete statistical systems would be encouraging. 

The main achievement of the present paper is to have shown that the Hamilton vector 
flow of CTM with arbitrary N gives a renormalization-group flow of random tensor networks, 
where the N = 2 case, in particular, corresponds to the Ising model on random networks. In 
the previous paper [3j, we considered the Hamiltonian of CTM, H = with “reasonable” 

choices of M a . Though it was successful in the N = 2 case, we have shown in this paper that 
the previous procedure of taking H does not work for general N, and have argued that the 
correct one is given by H = (pa/hLa-, where (p a are the integration variables for describing random 
tensor networks. Here an advantage of the present procedure is that H is uniquely determined 
by the dynamics of random tensor networks, but not by the ambiguous “reasonable” choices of 
the previous procedure. In fact, applied on random tensor networks, H = 0 a "H a is an operator 
which randomly inserts vertices on connecting lines, and therefore it increases sizes of tensor 
networks. This provides an intuitive understanding of the role of H as a renormalization 
process. We have performed the detailed analysis of the process, and have actually derived the 
Hamilton vector flow of CTM as a renormalization-group flow of the random tensor network. 
In the last section, we have proven a theorem which relates the phase transitions of the random 
tensor network and discontinuity of the renormalization-group flow. 

The renormalization-group flow which we have obtained has discontinuities on the first- 
order phase transition lines. However, there is a critical argument on whether a renormalization- 
group flow has discontinuities on a first-order phase transition line [31]. Since the argument 
basically assumes a regular lattice-like structure of a system, it would be interesting to inves¬ 
tigate a similar argument for a system with a random network structure. The random tensor 
network would give an interesting playground to deepen the idea of the renormalization group 
in a wider situation. 

Finally, let us comment on possible directions of further study, based on the achievements 
of the present paper. One is the classification of the fixed points of the Hamilton vector flow. 
This will provide the classification of the phases and their transitions of the random tensor 
network. This would also be interesting from the view point of quantum gravity. As discussed 
in [27], the physical wave functions of CTM may have peaks at the values of P invariant 
under some symmetries. In general, on such symmetric values of P, (f> may have multiple 
solutions, and therefore they are phase transition points. Such interplay between peaks and 
phase transitions may give interesting insights into quantum gravity. Another direction would 
be to pursue possible relations between the renormalization procedure of the random tensor 
network and that of the standard field theory. In fact, the “grand” partition function [3j of the 
random tensor network can be arranged to take the form of a partition function of field theory 
on a lattice by considering an index set labelling lattice points and taking P so as to respect 
locality. Then, the Hamilton vector flow of CTM may be regarded as a renormalization-group 
flow of the standard field theory. It would be highly interesting if CTM makes a bridge between 
quantum gravity and the standard field theory. 
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A Explicit expressions of the constraints 

In this appendix, we give the explicit expressions of the constraints, (TT2|) . (TT3|) . (fT4]h in the 
forms used in Section [4] for N = 2,3. 


A.l N = 2 


In a subspace, 


-Pm = 1, Pii2 = y , -P 122 = x \, -P 222 = x 2 , (64) 

with fixed y, which contains (1401) as a special case, the constraints are given by 


(fti,ft 2 ) 

= \ f 3 ^ 1 + y2 ^ + 3 ( x + x i)y^i~ + ( Xl + 2x2 + y( x 2 + 2 ^)) + y)-S~, 

6 V aiin ain 2 oxi ox 2 

3(1 + x 1 )y^~ + (xi + 2xl + y ( x 2 + 2y)) ^ 5 — + + y)jr~ + 3(x? + , 

Oriii dP 1U OX\ OX 2 ' 


Ji 2 OC -3 y 


d 


dP 


111 


+ (1-2x0 


d 


dP 


112 


V oc 


d 


d 


dP 


y 


in 


dP, 


X\ 


112 


d 

dxi 


d c) 

+ (-x 2 + 2 y)- -h 3xi — : 

OX i OX 2 

d 


X2 


dx 9 


(65) 

( 66 ) 

(67) 


A Hamilton vector flow is generated by (l20]h and, for a given Af a , the multipliers associated 
to the kinematical symmetries, J\f\ 2 and A7, can be determined so that the flow stays in the 
gauge AMD- Then, since P m and P\ 12 are kept constant by such a flow, determining such a 
Hamilton vector flow is actually equivalent to considering H = J\f a PL a , where Q p and Q p 
are substituted by solving the linear equations, Y12 = T> = 0. Here we do not write the explicit 
resultant expression of if, since it is rather long but the procedure itself is elementary. An 
important issue in the procedure is that there exist exceptional points characterized by 

3y 2 + 1 — 2xi = 0, (68) 
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where the set of linear equations, Y 12 = T> = 0 , become singular and cannot be solved for 
• qp , Q p • On these points, Af \2 and Af cannot be chosen so that the gauge be kept. These 
are the gauge singularities in Figj4j located at x\ = 0.635 for y = 0.3. 


A.2 N = 3 

The derivation of the Hamilton vector flow in the S 3 symmetric subspace (j4Tj) is similar to 
the N — 2 case. A difference is that, in (|20|) . the multiplier associated to J a b must be set 
Af a b = 0 to keep the S 3 invariance. For J\f a = 1, which is S 3 symmetric, one can choose A f 
appropriately to keep Pm = 1, and can obtain the Hamilton vector flow as 

H Ss = 2l 1 + X2 ^(1 + 3x! - 6xl + 2x ^^r + 6xi ( 1 _ ' ( 69 ) 

This is used to draw the Hamilton vector flow in Figl5] 


References 

[1] K. G. Wilson, “The Renormalization Group: Critical Phenomena and the Kondo Prob¬ 
lem,” Rev. Mod. Phys. 47 (1975) 773. 

[2] K. G. Wilson and J. B. Kogut, “The Renormalization group and the epsilon expansion,” 
Phys. Rept. 12 (1974) 75. 

[3] N. Sasakura and Y. Sato, “Ising model on random networks and the canonical tensor 
model,” PTEP 2014, no. 5, 053B03 (2014) [arXiv:1401.7806 [hep-th]]. 

[4] N. Sasakura and Y. Sato, “Exact Free Energies of Statistical Systems on Random Net¬ 
works,” SIGMA 10, 087 (2014) [arXiv:1402.0740 [hep-th]]. 

[5] N. Sasakura, “Canonical tensor models with local time,” Int. J. Mod. Phys. A 27 (2012) 
1250020 [arXiv:1111.2790 [hep-th]]. 

[ 6 ] N. Sasakura, “Uniqueness of canonical tensor model with local time,” Int. J. Mod. Phys. 
A 27 (2012) 1250096 [arXiv: 1203.0421 [hep-th]]. 

[7] N. Sasakura, “A canonical rank-three tensor model with a scaling constraint,” Int. J. 
Mod. Phys. A 28 (2013) 1 [arXiv:1302.1656 [hep-th]]. 

[ 8 ] J. Ambjorn, B. Durhuus and T. Jonsson, “Three-Dimensional Simplicial Quantum Grav¬ 
ity And Generalized Matrix Models,” Mod. Phys. Lett. A 6 , 1133 (1991). 

[9] N. Sasakura, “Tensor Model For Gravity And Orientability Of Manifold,” Mod. Phys. 
Lett. A 6 , 2613 (1991). 


21 







[10] N. Godfrey and M. Gross, “Simplicial Quantum Gravity In More Than Two-Dimensions,” 
Phys. Rev. D 43, 1749 (1991). 

[11] R. Gurau, “Colored Group Field Theory,” Commun. Math. Phys. 304 (2011) 69 
[arXiv:0907.2582 [hep-tli]]. 

[12] R. Gurau and J. P. Ryan, “Colored Tensor Models - a review,” SIGMA 8, 020 (2012) 
[arXiv: 1109.4812 [hep-th]]. 

[13] R. Gurau, “The 1/N Expansion of Tensor Models Beyond Perturbation Theory,” Com¬ 
mun. Math. Phys. 330, 973 (2014) [arXiv: 1304.2666 [math-ph]]. 

[14] V. Rivasseau, “The Tensor Track, III,” Fortsch. Phys. 62, 81 (2014) [arXiv:1311.1461 
[hep-th]]. 

[15] J. Ben Geloun, “Renormalizable Models in Rank d > 2 Tensorial Group Field Theory,” 
Commun. Math. Phys. 332, 117 (2014) [arXiv:1306.1201 [hep-th]]. 

[16] D. Oriti, “Group Field Theory and Loop Quantum Gravity,” arXiv: 1408.7112 [gr-qc]. 

[17] D. Ousmane Samary, “Beta functions ofU(l) d gauge invariant just renormalizable tensor 
models,” Phys. Rev. D 88, no. 10, 105003 (2013) [arXiv:1303.7256 [hep-th]]. 

[18] A. Eichhorn and T. Koslowski, “Continuum limit in matrix models for quantum grav¬ 
ity from the Functional Renormalization Group,” Phys. Rev. D 88, 084016 (2013) 
[arXiv:1309.1690 [gr-qc]]. 

[19] M. Raasakka and A. Tanasa, “Next-to-leading order in the large N expansion of the 
nmlti-orientable random tensor model,” arXiv: 1310.3132 [hep-th]. 

[20] S. Carrozza, “Discrete Renormalization Group for SU(2) Tensorial Group Field Theory,” 
arXiv:1407.4615 [hep-th], 

[21] J. Ben Geloun and R. Toriumi, “Parametric Representation of Rank d Tensorial Group 
Field Theory: Abelian Models with Kinetic Term |Ps| +/b” arXiv:1409.0398 [hep-th], 

[22] V. A. Nguyen, S. Dartois and B. Eynard, “An analysis of the intermediate field theory 
of T 4 tensor model,” arXiv: 1409.5751 [math-ph], 

[23] D. Benedetti, J. Ben Geloun and D. Oriti, “Functional Renormalisation Group Approach 
for Tensorial Group Field Theory: a Rank-3 Model,” arXiv:1411.3180 [hep-th], 

[24] T. Delepouve and V. Rivasseau, “Constructive Tensor Field Theory: The Tg Model,” 
arXiv:1412.5091 [math-ph], 

[25] N. Sasakura and Y. Sato, “Interpreting canonical tensor model in minisuperspace,” Phys. 
Lett. B 732, 32 (2014) [arXiv:1401.2062 [hep-th]]. 

[26] N. Sasakura, “Quantum canonical tensor model and an exact wave function,” Int. J. Mod. 
Phys. A 28 (2013) 1350111 [arXiv:1305.6389 [hep-th]]. 


22 


[27] G. Narain, N. Sasakura and Y. Sato, “Physical states in the canonical tensor model from 
the perspective of random tensor networks,” JHEP 1501, 010 (2015) [arXiv:1410.2683 
[hep-tli]]. 

[28] See, S. N. Dorogovtsev, A. V. Goltsev, J. F. F. Mendes, “Critical phenomena in complex 
networks,” Rev. Mod. Phys. 80, 1275 (2008) [arXiv:0705.0010 [cond-mat]], and references 
therein for the dynamics of statistical systems on random networks. 

[29] R. L. Arnowitt, S. Deser and C. W. Misner, “Canonical variables for general relativity,” 
Phys. Rev. 117, 1595 (1960); R. L. Arnowitt, S. Deser and C. W. Misner, “The Dynamics 
of general relativity,” arXiv:gr-qc/0405109, 

[30] For instance, see N. Goldenfeld, “Lectures on Phase Transitions and the Renormalization 
Group,” ISBN-020155408955408. 

[31] For instance, see A. D. Sokal, A. C. D. van Enter and R. Fernandez, “Renormalization 
transformations in the vicinity of first-order phase transitions: What can and cannot go 
wrong,” Phys. Rev. Lett. 66, 3253 (1991), and references therein. 


23 


